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ABSTRACT 

We have carried out a submillimeter continuum and spectroscopic study of the W43 main complex, a 
massive star-forming region, which harbors a giant H n region. The maps reveal a filamentary structure 
containing ~ 50 fragments with masses of 40 — 4 000 M and typical diameters of 0.25 pc. Their large 
sizes, large non-thermal velocities (Av ~ 5 kms" 1 ), and high densities (n H2 ~ 10 6 cm -3 ) suggest that 
they are protoclusters and excellent sites to form massive stars. Follow-up observations are necessary, 
but we have already identified three protoclusters to be very good candidates for containing very young 
massive protostars. The starburst cluster, that excites the giant H n region has a large impact on the 
molecular complex. However, it remains unclear if this first episode of star formation is triggering the 
formation of new massive stars, through ionization shocks crossing the closeby molecular clouds. W43 
is thus an ideal laboratory to investigate massive star formation from the protostellar phase to that of 
giant H II regions. Moreover, the very active star-forming complex W43 may be considered a Galactic 
mini-starburst region that could be used as a miniature model of starburst galaxies. 

Subject headings: dust — H n regions 
— submillimeter 



ISM: individual (W43) — ISM: structure — stars: formation 



1. INTRODUCTION 



High-mass (OB; M* > 8 M Q ) stars are believed to form 
in clusters within molecular cloud complexes. Due to their 
high ultraviolet luminosity, massive young stellar objects 
(YSOs) first heat, then ionize, and disrupt their surround- 
ing molecular cloud. Owing to their large distances to 
the Sun and the complex interplay between massive YSOs 
and the neighboring interstellar medium, the formation of 
high-mass stars is still poorly understood. Many studies 
have been devoted to embedded YSOs that have already 
developed an H n region and thus are easily detectable 
in far-infrared and centimeter continuum surveys (see a 
review by Churchwell 1999). In contrast, possible precur- 
sors of ultracompact H n regions (UCH us), i.e. massive 
YSOs in their main building phase, have only been dis- 
covered recently (Hunter et al. 2000; Brand et al. 2001; 
Sridharan et al. 2002 and references therein). These mas- 
sive protostars, also called high-mass protostellar objects 
(HMPOs), are inconspicuous in the IRAS bands and have 
no, or very little free-free emission. They can be detected 
in (sub)millimeter dust continuum and molecular high- 
density tracers and frequently display H2O/CH3OH maser 
emission. Therefore, submillimeter continuum imaging of 
molecular cloud complexes is ideal to probe clouds sur- 
rounding UCH us and HMPOs simultaneously. Making a 
census of those deeply embedded phases is the first neces- 



sary step to gain insight into the processes leading to the 
formation of a massive star. 

The W43 star-forming complex is located in the inner 
spiral arm of our Galaxy, at 5.5 kpc from the Sun (Wilson 
et al. 1970). W43 is well-known for its giant H 11 region 
emitting 10 51 Lyman continuum photons per second and 
a far-infrared continuum luminosity of ~ 3.5 x 10 6 Lq 
(Smith, Biermann, & Mezger 1978; Lester et al. 1985). 
The main ionizing source of W43 was discovered in near- 
infrared images by Lester et al. (1985) and confirmed by 
Blum, Damineli, & Conti (1999) as a cluster of Wolf-Rayet 
(WR) and OB main sequence stars. The inner 20 pc of 

W43 (at I ~ 30.75°, b 0.06°) contain this WR/OB 

cluster and a 10 6 M Q molecular cloud called G30. 8-0.0. 
The W43 main complex was mapped in CO lines and parts 
of it in higher density tracers such as H 2 CO and HCO+ 
lines, and 1.3 mm continuum (Bieging, Wilson, & Downes 
1982; Liszt 1995; Mooney et al. 1995). Several sources were 
identified by these authors but a complete and comprehen- 
sive sample of UCH n regions and massive protostars is 
still lacking. Identifying such a sample is essential for the 
present paper, which aims at presenting a global scenario 
of (massive) star formation in W43. Therefore, we will 
carefully investigate the distribution of dense molecular 
clouds in the W43 main complex, in particular to identify 
sites of future or on-going massive star formation. 
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A detailed study of W43 may help constrain the prop- 
erties of distant starburst galaxies. Indeed, the large lumi- 
nosity and ionizing flux of the W43 giant H n region are 
similar to those of NGC 3603 or M17 (10 5 - I0 7 L Q and 
10 50 — 10 51 Lyes -1 ) , taken to be representative of clus- 
ters and H II regions in starburst galaxies (e.g. Tapia et al. 
2001). The NGC 3603 stellar cluster is qualified as "star- 
burst" because it consists of several tens of WR, O, and 
B-type stars (e.g. Brandl et al. 1999). The W43 cluster of 
main sequence stars is likely to be similar but more heav- 
ily reddened (A v = 30 mag versus 4 mag in NGC 3603). 
Since the W43 stellar cluster is closely associated with gi- 
ant molecular clouds, we may be witnessing another burst 
of star formation. This would make the W43 main complex 
a Galactic mini-starburst region, i.e. a miniature model of 
the stellar and gas content of starburst regions in distant 
galaxies. 

In the present paper, we report a submillimcter contin- 
uum and spectroscopic study of the W43 main complex. 
From the imaging and deep spectroscopic measurements 
presented in Sect. 2, we make a complete census of the 
compact and dense cloud fragments in W43 (Sect. 3). In 
Sect. 4, we investigate the impact of the giant H n region 
on the molecular complex, the nature of these cloud frag- 
ments, and the global characteristics of the mini-starburst 
W43. Finally, Sect. 5 summarizes our conclusions. 

2. OBSERVATIONS 

2.1. Dust Continuum Observations 

We mapped the dust continuum emission of the W43 
main complex at 1.3 mm and 350 fim in February and 
April 1999. We used the 37-channel MPIfR MAMBO 
bolometer array installed at the IRAM 2 30-meter telescope 
on Pico Veleta (Spain) to make the 1.3 mm map. We per- 
formed the 350 /im observations at the CSO 3 10.4-meter 
radiotelescope on Mauna Kea (Hawaii) , equipped with the 
24-pixel SHARC bolometer camera. The passband of the 
MAMBO (respectively SHARC) bolometer array has an 
equivalent width of « 70 GHz (resp. 100 GHz) and is cen- 
tered at v cS « 240 GHz (resp. 860 GHz) (Kreysa et al. 
1998; Hunter, Benford, & Serabyn 1996). 

Both 1.3 mm and 350 /im images cover ~ 12' x 10' and 
were obtained by combining partially overlapping on-the- 
fly maps. In the dual-beam on-the-fly mapping mode, the 
telescope is scanned continuously in azimuth along each 
row, while wobbling. For each channel, the raw data cor- 
responding to a single on-the-fly coverage consist of several 
rows taken at a series of elevation offsets. For MAMBO 
(respectively SHARC) maps, we used a scanning velocity 
of 4" sec -1 (resp. 10" sec -1 ), and a sampling of 2" (resp. 
4") in azimuth and 4" in elevation. The wobbler frequency 
was set to 2 Hz (resp. 4 Hz) and the wobbler throw in az- 
imuth was 70" (resp. 135"). The typical azimuthal size of 
individual maps was ~ 9' with MAMBO versus ~ 12' with 
SHARC. The MAMBO maps were reduced with the IRAM 
software for bolometer- array data (NIC; cf. Broguiere, 
Neri, & Sievers 1995) while the SHARC data were pro- 
cessed with the standard CSO programs CAMERA and 
REGRID. The IRAM and CSO programs both use the 
EKH restoration algorithm (Emerson, Klein, & Haslam 

2 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 

3 The CSO is operated by the California Institute of Technology under 



1979). 

The size of the main beam was measured to be 
HPBW ~ 11" for MAMBO maps and - 11.5" for SHARC 
maps using Uranus and Mars. The absolute pointing of 
each telescope was found to be accurate to within ~ 5". 
The data were taken with moderate weather conditions. 
The zenith atmospheric optical depth at 240 GHz (re- 
spectively 225 GHz) varied between ~ 0.2 and ~ 0.4 for 
MAMBO observations and was roughly ~ 0.06 for SHARC 
observations. Uranus and Mars were used for flux calibra- 
tion and the overall, absolute calibration uncertainty is 
estimated to be ~ 20% for both MAMBO and SHARC 
maps. 

2.2. Molecular Line Observations 

In June 2000 and 2001, we mapped the HCO+(3-2) ro- 
tational line emission of the W43 main complex and con- 
ducted deep integrations at some of the most prominent 
dust maxima in the (3-2) transitions of HCO + and its 
H 13 CO+ isotopomer. We used the facility 230 GHz re- 
ceiver of the CSO and its 50 MHz bandwidth acousto- 
optical spectrometer with a ~ 0.16 kms -1 spectral resolu- 
tion. The CSO FWHM angular resolution for frequencies 
of 267.5576 GHz and 260.2555 GHz is 28" - 29". 

We made two overlapping on-the-fly maps that cover 
a 8' x 10' area. Each of these maps consists of several 
rows continuously scanned in right ascension and taken at 
a series of declination offsets. The scanning velocity was 
7" sec -1 and the sampling 10" in both directions. The 
off position was taken at each declination offset, 20' away 
in right ascension. Besides this map, we performed deep 
measurements in position switching mode, using off posi- 
tions ±20' in azimuth. Pointing was checked on planets 
when available or W-AQL otherwise and was found to be 
stable to within 5" - 10". 

The data were reduced with the IRAM software for spec- 
tral lines (CLASS), applying low order baselines and av- 
eraging individual spectra with rms weighting. Our line 
observations were corrected for atmospheric attenuation, 
ohmic losses, rearward spillover and scattering, using the 
standard chopping wheel technique. The 225 GHz zenith 
optical depth was ~ 0.075 during the maps and varying 
between 0.08 and 0.18 during the position switch mea- 
surements. To convert the resulting antenna tempera- 
tures (T*) to main beam brightness temperatures (T MB ), 
we have divided the data by the main beam efficiency de- 
termined from observations of Mars, ?7 MB = 0.65. The 
calibration accuracy is ~ 25%. 

3. RESULTS AND ANALYSIS 

Our 1.3 mm and 350 /xm continuum maps are presented 
in Figs. la-b. They cover a 17 pc x 20 pc area, roughly 
centered on the stellar WR/OB cluster. The HCO+(3-2) 
map covers a slightly smaller area (14 pc x 16 pc) and 
is shown in Figs. 2a- f, where the integrated intensity is 
displayed along with five velocity channels ranging from 
86 kms -1 to 102 kins" 1 . Figs, la-b and 2a all show that 
the dense molecular gas in the W43 main complex follows 
a global "Z" -shaped filamentary structure. Some more 
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compact fragments are observed within this filament. We 
make a complete census of these fragments in Sect. 3.1 and 
estimate their temperature, mass, and other characteris- 
tics in Sect. 3.2. 

3.1. Census of Compact Cloud Fragments in W43 
3.1.1. Cloud structure 

Our immediate goal is to investigate the submillimeter 
continuum maps of W43 to identify cloud fragments which 
could form massive stars. The maps shown in Figs, la-b 
are in very good agreement with each other (except for the 
north-western corner, see below) and mainly trace the op- 
tically thin thermal emission of dust in molecular clouds. 
Examination of Figs, la-b reveals a wealth of cloud struc- 
tures with size scales ranging from 11" (~ 0.3 pc) to more 
than 500" (~ 13 pc). The dynamical range within these 
maps is also very high: 450(7 in the MAMBO map and 
75a in the SHARC map. To assess the complexity of 
the submillimeter emission of W43, we use the Gauss- 
clumps program (Kramer et al. 1998), which we have 
adapted to continuum observations. The emission above 
5cr = 60 Jybeam -1 in Fig. la can fairly well be described 
by a sum of ~ 300 cloud fragments with 2-D Gaussian 
shapes. Since most of such fragments consist of moderate- 
density gas (n H2 < 10 4 cm -3 , see below), they are unlikely 
sites for massive star formation. 

We thus need to separate compact and dense cloud frag- 
ments from their surrounding, lower density clouds. We 
use the method developed by Motte & Andre (2001a) to 
estimate the outer diameters of the four main submillime- 
ter fragments (W43-MM1 to MM4 labelled in Fig. la). 
In both MAMBO and SHARC maps, we measure outer 
diameters of ~ 40" corresponding to ~ 1 pc. We do 
not expect the weaker compact fragments to have the ex- 
act same radial extent but, for simplicity, we use a 40" 
lengthscale to extract all dense fragments. We stress that 
such a cloud decomposition is not unique as these ~ 1 pc 
cloud fragments will probably fragment further to form 
more than one single star. It permits one to select struc- 
tures with peak density larger than ~ 10 cm~ 3 , when a 
5(7 = 60 mjy/ll"-beam 1.3 mm flux, a 20 K dust temper- 
ature and a k 13 mm = 0.01 cm 2 g _1 dust mass opacity are 
assumed (calculated from e.g. Eq. (1') of Motte, Andre, 
& Neri 1998). 

We measure the submillimeter flux of each compact frag- 
ment in two steps. We first filter out all spatial scales 
larger than 40" from the MAMBO map and identify 51 
fragments above the 5cr level. We then derive the 1.3 mm 
Gaussian characteristics of each fragment by applying an 
improved version of the technique we developed for our 
previous studies (Motte et al. 1998, 2001). It uses a mul- 
tiresolution analysis to refine the extraction of the compact 
fragment from its local background emission (see Starck, 
Murtagh, & Bijaoui 1998; Motte & Andre 2001b) and 
the Gaussclumps program to fit its Gaussian parameters 
(Kramer et al. 1998). The fitted fragments are indicated 
as black ellipses in Fig. la and account for ~ 15% of the 
total 1.3 mm flux measured. We perform the same analysis 
on the SHARC map, where 23 of the 51 1.3 mm fragments 
are recovered above 3(7 = 9 Jybeam -1 and within 10" of 
the 1.3 mm positions. 

These compact submillimeter fragments have basic char- 



acteristics listed in Table 1. Col. 1 gives their adopted 
names (where a small number corresponds to a high level of 
detection during the first step of our analysis) and Cols. 2-3 
their coordinates. The Gaussian peak fluxes, FWHM sizes, 
and integrated fluxes measured at 1.3 mm are listed in 
Cols. 4-6, along with the integrated fluxes of their 350 /xm 
counterparts in Col. 7. The spectral indices a^°° (where 
S v oc v a ) we measure from the integrated fluxes S^f mm 
and S 3 f* are given in Col. 8. The error bar on aH°° 
varies from ±0.1 up to ±0.4 for the weaker fragments. 

The next step of the analysis is to ensure that the sub- 
millimeter fragments selected above are indeed dense cloud 
fragments. We thus need to check the contamination of 
submillimeter fluxes by lines and free-free emission. Con- 
tamination by CO line emission is negligible for those frag- 
ments which are dense and compact. In contrast, the "hot 
core" line emission displayed by several fragments (Motte 
et al. in prep.) contributes to the submillimeter contin- 
uum fluxes we measure. Such "hot cores" are pointlike (in 
a 0.3 pc beam) and embedded in dense fragments, which 
have a S^JS*^ ratio of - 2 (cf. Table 1). Their 
emission would thus very unlikely dominate the submil- 
limeter continuum. We consider that 25% is an upper 
limit for the line contamination of the W43 submillimeter 
fragments. 

Without a doubt, the main source of contamination of 
1.3 mm fluxes is the free-free emission from H n regions. 
To estimate such a contamination, we compare our 1.3 mm 
and 350 /im images with centimeter free-free data from 
the literature. Strikingly, the maps of Figs, la-b differ 
in their north-western corner where the 1.3 mm arc-like 
structure disapears at 350 /im. For a further comparison 
we compute the spectral index map between 350 /im and 
1.3 mm (see Fig. 3). Most of the map shows an index 
of 

a 35o° ~ 4 ± 0.5, suggestive of thermal emission of dust 
at > 30 K. In contrast, at the location of the giant H II 
region (Lester et al. 1985; Liszt, Braun, & Greisen 1993; 
Balser, Goss, & De Pree 2001), the map has a lower index: 
aH°° <~ 2.5 — 3.5. Since there is a good correspondance be- 
tween these low aH°° values and the large-scale ionized 
gas (cf. Fig. 3), the 1.3 mm fluxes measured in that re- 
gion are likely contaminated by free-free emission. Most 
of the free-free emission in W43 arises from its giant H II 
region and is probably optically thin at centimeter wave- 
lengths. Any ultra-compact H II regions in W43 should 
have a free-free emission at least partly optically thick. 
The spectral indices computed for the two most compact 
centimeter sources are however indicating that their free- 
free emission is mostly optically thin when A < 3 cm (see 
Sect. 3.2.4). We thus hereafter assume that all the 3.5 cm 
emission mapped in W43 at ~ 10" angular resolution is 
optically thin. 

We estimate the contribution of free-free emission to the 
1.3 mm fluxes measured in Table 1 by using the 3.5 cm 
map of Balser et al. (2001), which has a similar angular 
resolution to our continuum data (9" vs. 11"). The 21 cm 
map of Liszt et al. (1993) is used when the fragments are 
lying outside the 3.5 cm primary beam. We filter out all 
spatial scales larger than 40" and integrate the centimeter 
emission at the position of each 1.3 mm fragment. Then 
we assume a S v oc v~ 0A spectra typical for optically thin 
free-free emission to estimate the flux of the ionized gas 
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at 1.3 mm (cf. Col. 9 in Table 1). Within the giant H n 
region powered by the WR/OB cluster, the free-free con- 
tribution to the twenty 1.3 mm compact fragments (in- 
cluding W43-MM6, MM8, MM14, MM11, MM4, MM20, 
MM15 and MM25 shown in Fig. la) goes up to 70% with 
an average value of 16%. Outside this region, the free-free 
contribution is lower than 10%, with the notable exception 
of W43-MM13, which coincides with a centimeter source 
(cf. Fig. 3). Altogether the level of free- free contamination 
is low for the compact continuum fragments identified in 
Fig. la. 

In the following, we subtract the free-free emission esti- 
mated above from the 1.3 mm fluxes of Table 1. Three of 
the 51 fragments (namely W43-MM32, MM46, MM50) are 
falling below the 5c limit we set for their initial detection 
and are removed from the sample. We thus end up with 48 
1.3 mm fragments, 22 of which are also detected at 350 /zm. 
We are confident that the 48 compact fragments identified 
in Table 3 represent genuine dust continuum sources. 

3.1.2. Cloud kinematics 

We make here a first order kinematics analysis of the 
dense cloud fragments observed in the W43 main molecu- 
lar complex. The HCO + (3-2) emission in Fig. 2a displays 
the same filamentary structure as the dust continuum 
maps in Figs. la-b. We selected several compact dust frag- 
ments at various locations in the W43 main complex and 
made deep integrations in HCO+(3-2) and H 13 CO+(3-2). 
In Table 2, for each submillimcter fragment (Col. 1) and 
optically thin molecular transition (Col. 2), we separate 
the line into one to four velocity components (see Col. 3). 
We choose this decomposition to better describe the lines 
which are broad and structured but none of those compo- 
nents are perfectly fitted by a Gaussian. We give in Cols. 4- 
7 of Table 2 the Gaussian parameters of the fitted com- 
ponents, i.e. local velocity at rest, peak line temperature, 
integrated intensity and line width. Figs. 4a-e display the 
spectra taken at some of these locations. HCO + emission 
has been successfully detected for each selected submil- 
limcter fragment. Even 1.3 mm locations W43-MM15&20 
and MM25, where heating from the WR/OB cluster and 
contamination by free-free emission is highly probable, 
contain high-density gas. As suggested by a careful exam- 
ination of Fig. 2a , the HCO+(3-2) line emission of W43- 
MM1, MM2, MM3, MM9, and MM10 is optically thick 
(cf. Figs. 4a and 4c). The HCO+(3-2) emission of the 
other submillimcter fragments is probably optically thin, 
but their lines often display a complex shape: from mul- 
tiple velocity components (W43-MM15&20, MM4, MM7, 
cf. Figs. 4c-d) to a single component with a broad non- 
Gaussian line shape (W43-MM6&8, MM25, cf. Fig. 4b). 
We stress that the more complex is the line, the closer in 
projection to the WR/OB cluster is the fragment. De- 
spite those limitations, we will use the values of v LSR and 
Av given in Tabic 2 to further analyze the molecular cloud 
kinematics. 

As illustrated in Figs. 2b-f and Figs. 4a-d, there is 
a velocity gradient from the southern part of the fila- 
ment (f LS R ~ 85 kms -1 ) to its northern part (f LSR ~ 
105 kms -1 ). Liszt (1995) observed similar gradients in 
HCO+(1-0) and 13 CO(1-0) and proposed that the south- 
ern gas is in front of W43 and the northern gas behind (see 



also Bieging et al. 1982). This velocity gradient roughly 
develops along the Galactic plane and encompasses the re- 
combination line velocity found for the ionized gas in front 
of the WR/OB cluster (92.4 kms -1 according to Balser et 
al. 2001). It thus suggests that the starburst cluster is inti- 
mately associated with the W43 main molecular complex. 

3.2. Characteristics of the Compact Cloud Fragments 

We derive the basic properties of the compact (< 1 pc) 
cloud fragments of W43 using their dust continuum and 
HCO + line parameters given in Tables 1-2. Table 3 lists 
the fragment name (Col. 1), the adopted dust tempera- 
ture (Col. 2), the estimated submillimcter mass (Col. 3), 
gas density (Col. 4), and virial mass (Col. 5). In the fol- 
lowing, we describe the methods and assumptions used 
to estimate those characteristics and we search for stellar 
activity signatures coincident with the dense cloud frag- 
ments. 

3.2.1. Dust temperature 

The temperature is unlikely to be homogeneous in the 
W43 molecular complex because it is known to form mas- 
sive stars. As a matter of fact, Lester et al. (1985) mea- 
sured dust temperatures ranging from 40 K to 90 K, using 
50 /zm and 100 /im images made at the Kuiper Airborne 
Observatory (KAO). We expect the temperature of W43 
compact fragments to be lower as they are largely em- 
bedded within clouds. The first confirmation comes from 
the spectral indices of the 22 fragments detected at both 
1.3 mm and 350 ^tm. When they are corrected for free- 
free contamination, these indices have a median value of 
a 1 ™ - 33 ( see Cols - 8-9 of Table 1). For a pure ther- 
mal and optically thin dust emission with a standard dust 
opacity index of = 2 (where k v cx v 13 ), this corresponds 
to the averaged dust temperature of <~ 20 K. 

In an effort to better constrain the average temperature 
of the W43 compact fragments we compare their spec- 
tral energy distribution (SED) to graybody models (cf. 
Figs. 5a-b). Their SED is built from the 1.3 mm and 
350 /im integrated fluxes measured in this paper (cf. Ta- 
ble 1) and the 100 /im and 50 /im fluxes read from Fig. 1 of 
Lester et al. (1985) for pointlike sources (in a 50" beam). 
The fluxes of MSX sources located close to the submil- 
limcter fragments or corresponding upper limits are dis- 
played but not used (see Sect. 3.2.4). Our SED analysis 
shows that graybody models with dust opacity index of 
= 2 might be more appropriate than those with = 1.5 
(where n v oc v^, cf. Fig. 5a). This result agrees with 
the multiwavelength submillimeter studies of other mas- 
sive YSOs (e.g. Hobson et al. 1993). We therefore use 
graybody models with — 2, while knowing that temper- 
atures increase by only a few degrees when = 1.5 (e.g. 
Fig. 5a). Since mid- to far- infrared fluxes are often missing 
and upper limits remain high (see also Sect. 4.1.1), we only 
successfully constrain the dust temperature of four dense 
fragments: W43-MM1 (T dust - 19 K), MM2 (~ 23 K), 
MM3 (< 20 K), and MM7 (< 22 K). We also derive the 
bolometric luminosity and the submillimeter to bolomet- 
ric luminosity ratio of W43-MM1 (L bol ~ 2.3 x 10 4 L , 
L x > 350 ^/L bol ~ 3%), MM2 (~ 2.4 x 10 4 L Q , ~ 1.5%), 
MM3 (< 1.2xl0 4 L Q , > 2.5%), and MM7 (< l.OxlO 4 L , 
> 1.5%). Despite its incompleteness, the present analysis 
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suggests that most compact cloud fragments in W43 are 
rather cold (T dust ~ 20 K) and have high bolometric lu- 
minosities (L bol <~ 10 4 Lq). Their masses are estimated 
below. 

3.2.2. Mass estimated from the submillimeter continuum 

We believe that the submillimeter continuum emission 
of cloud fragments identified in Sect. 3.1.1 is mainly ther- 
mal dust emission, which is largely optically thin. For any 
given dust properties and gas-to-dust ratio, the 1.3 mm 
and 350 /im fluxes are thus directly related to the total 
(gas + dust) mass of the fragments. For present mass es- 
timates, we use the 1.3 mm integrated fluxes (5 , 1 i ° t mm , see 
Table 1) because they provide more homogeneous mea- 
sures. We correct those fluxes for any residual free-free 
contamination (cf. Col. 9 in Table 1) and derive the "sub- 
millimeter mass" (M Bmm , given in Table 3), as follows: 

j^j- _ (S 13 mm ) d 

1^1.3 mm B 1 3 mm (T du!t ) 

X (o.Olcm 2 g- 1 ) (20 k) ' ^ 

where k 13 mm is the dust opacity per unit mass column 
density at 1.3 mm, and B 13 mm (T duat ) is the Planck func- 
tion for a dust temperature T dUBt . 

The dust mass opacity (including dust properties and 
gas-to-dust mass ratio) is likely to vary with density, 
temperature, and the evolutionary state of the emitting 
medium (Hcnning, Michel, & Stognicnko 1995). Models of 
dust in low-mass protostellar cores (e.g. Ossenkopf & Hcn- 
ning 1994) suggest that a value of k 13 mm = 0.01 cm 2 g -1 
is well suited for cool (10 — 30 K) and high-density 
(n H 2 10 5 cm -3 ) cloud fragments. This value agrees 
with the recent cross-comparisons of dust emission sur- 
rounding an UCH 11 region with its CO ice absorption 
and gas emission (van der Tak 2002). We thus choose a 
dust opacity per unit (gas + dust) mass column density 
of Ki.3mm = 0.01 cm 2 g _1 . We use this value for all the 
dense fragments because the average dust properties of 
these ~ 0.25 pc structures should not change drastically 
when they contain a massive protostar or an UCH II. We 
estimate that the absolute value taken for the dust mass 
opacity is uncertain by at least a factor of 2. 

The temperature to be used in Eq. (1) is the mass- 
weighted dust temperature of the cloud fragments, its 
value can be determined from the gray-body fitting of their 
SED. In Sect. 3.2.1, we recommend 20 K for most of the 
dense cloud fragments. We thus assume T dust = 20 K for 
the fragments deeply embedded into the "Z" -shaped fila- 
ment and T dust = 30 K for those that may be less shielded 
against the radiation of the WR/OB cluster (see Col. 2 of 
Table 3). 

3.2.3. Virial mass 

We estimate the virial mass of the W43 compact cloud 
fragments assuming that they are spherical and have a 

4 See http://www.ipac.caltech.edu/ipac/msx/msx.html. 

5 http://www.ipac.caltech.edu/ipac/iras/hires_over.html 



p oc r 2 density distribution: M vir = 3Ra 2 /G (e.g. 
Bertoldi & McKee 1992, cf. Table 3). We use the line 
width (Aw) listed in Table 2 and the diameter measured 
at 1.3 mm (FWHM in Table 1). We subtract the ther- 
mal component of the line width using the temperature 
adopted in Table 3. As the thermal component of the 
velocity dispersion is rather small (a e ~ 0.3 kms -1 when 
Av ~ 5 kms -1 ), the error introduced by assuming that the 
kinetic temperature equals the dust temperature is negli- 
gible. Whenever the optically thin line consists of several 
velocity components, we use their median line width and 
indicate the number of components in Col. 5 of Table 3. 

The targetted positions for HCO + measurements gener- 
ally correspond to a single submillimeter fragment embed- 
ded in a lower-density filament (cf. Sect. 3.1.1). Therefore 
the line width and then the virial mass derived from these 
measurements are upper limits for the width and mass of 
the submillimeter fragments alone. Follow-up maps with 
higher density tracers are needed to select the gas only 
associated with the submillimeter fragment and estimate 
a more accurate virial mass. 

In a few cases, HCO + lines are optically thin but re- 
veal several velocity components. The comparison with 
the line spectra of the H 13 CO + isotopomer can help iso- 
late the velocity component most probably associated with 
the compact and dense dust fragment. This is the case of 
W43-MM4 and MM7 where the low velocity components 
(at w LSR ~ 91.8 kms -1 and ~ 90.6 kms -1 ) have larger 
H 13 CO + /HCO + intensity ratios and should thus consist 
of higher density gas (cf. Table 2 and Fig. 4d). When we 
cannot state that a submillimeter fragment corresponds 
to a single velocity component, we simply assume that it 
consists of several cloud fragments with the same mass. 

According to Pound & Blitz (1993), a fragment remains 
gravitationally bound as long as M/M vir > 0.5. Since the 
submillimeter masses are uncertain by a factor of 2 and 
the virial masses are upper limits, we estimate that the 
M Bmm /M vir ratio must be at least larger than 0.2 to sug- 
gest gravitational boundedness. 

3.2.4. Coincidence with signposts of stellar activity 

Various published catalogues can be used to search for 
stellar activity within the dense cloud fragments listed in 
Table 3. In addition to the IRAS survey, the Midcourse 
Space Experiment (MSX 4 ) provides pointlike sources ob- 
served at 8, 12, 15, and 21 /im with a 20" resolution 
(Egan et al. 2001). MSX G30.7196-0.0854 and G30.6882- 
0.0726 (or IRAS 18450-0205) are the only two infrared 
sources relevant to our study because they coincide with 
high-density clouds and have < 40" sizes (see Fig. 6 and 
Sect. 4.1.1). MSX G30.7196-0.0854 may be associated 
with W43-MM3, but its nominal position is further away 
than what the maximum pointing uncertainties allow. As 
for MSX G30.6882-0.0726, it coincides with W43-MM13 
well within the 9" maximum positional errors. The 60 /im 
and 100 /im fluxes of the associated IRAS source are con- 
fused by the far-infrared emission of the W43 giant H 11 
region. With HIRES 5 (High Resolution Processing using 
the Maximum Correlation Method) images, we success- 
fully measure a 60 /im flux of 1 000 ± 500 Jy within a 
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~ 60" beam. 

Besides infrared sources, some small-diameter (< 40") 
radio sources are identified at 21, 6 and 3.5 cm by Zoone- 
matkermani et al. (1990), Garwood et al. (1988), Becker 
et al. (1994), and Balser et al. (2001). They roughly co- 
incide with W43-MM3, MM4, MM13, MM14 and MM20 
(see Fig. 3). W43-MM3 and MM4 are clearly the most 
compact of these sources with diameters of ~ 3" — 5". 
Their spectral indices measured from 21 to 6 cm and then 
from 6 to 3.5 cm are 0.5 and < 0.1 for W43-MM3 and 
~ 0.6 and < 0.5 for W43-MM4. Their free-free emission 
is thus partially optically thick from 21 to 3.5 cm and be- 
comes mostly optically thin at wavelengths shorter than 
3.5 cm. 

Furthermore, complete surveys of methanol, water and 
hydroxyl maser emission have been performed in W43 
(Walsh et al. 1998; Valdettaro et al. 2001; Braz & Epchtein 
1983). The CH 3 OH survey made at 6.7 GHz with ATCA 
has position uncertainties smaller than 2" (Walsh et al. 
1998) and the H2O maser positions observed at 22 GHz 
with the Medicina telescope are only secure to 25" (Valdet- 
taro et al. 2001). The OH type-I masers detected at 1.6- 
1.7 GHz and compiled by Braz & Epchtein (1983) have 
a positional accuracy of 5" or 30" (south-west of W43- 
MM2). In Fig. lb, W43-MM1 and MM11 are coincident 
with the CH3OH, H 2 and OH maser sources called Main3 
and Mainl, respectively. W43-MM2 is associated with 
CH3OH and OH masers while W43-MM6 and MM7 coin- 
cide with two H 2 maser sources, the second being called 
Main2. The position and systemic velocity of the methanol 
masers are well confined within the radial and velocity ex- 
tent of the compact fragments (< 10" and < 3 kms -1 
from their peaks, Walsh et al. 1998; see also Caswell et 
al. 1995). The water maser coincident with W43-MM1 
peaks 0.7 kms -1 away from its systemic velocity. That 
associated with W43-MM11 has one out of two velocities 
recorded within 2 kms -1 of the molecular emission peak 
(cf. Valdettaro et al. 2001 and references therein). These 
infrared, centimeter and maser sources are stellar activity 
signposts which will help investigate the stellar content of 
cloud fragments listed in Table 3. 

4. DISCUSSION 

We have made a complete census of cloud fragments 
with density greater than 10 4 cm -3 and size smaller than 
1 pc in W43. Before we discuss their characteristics and 
evolutionary state (in Sect. 4.2), we investigate the likely 
effects of the WR/OB cluster on the W43 molecular cloud 
complex (in Sect. 4.1). 

4.1. The Starburst Cluster Interacting with W43 
Molecular Clouds 

According to the near-infrared spectroscopic study of 
Blum et al. (1999), W43 harbors a cluster of main sequence 
stars that contains at least one Wolf-Rayet star and two 
O-type (super)giants. This young cluster is substantially 
more luminous (~ 3.5 x 10 6 Lq) than a typical OB associ- 
ation and qualifies as a starburst cluster. We first describe 
the characteristics of the giant H 11 region excited by this 
starburst cluster and then discuss its impact on the W43 
molecular clouds. 



4.1.1. The giant Hn region 

Previous centimeter continuum studies showed that the 
far-ultraviolet photons of the WR and OB stars in the 
W43 cluster created a large bubble of free-free emission 
also called giant H 11 region (see Balser et al. 2001 and 
references therein). The white dashed ellipse plotted in 
Fig. 3 outlines this bubble of ionized gas, which has a full 
size of ~ 4 pc x 6 pc. The north-eastern and south-western 
extensions of the W43 free-free emission may also be as- 
sociated with gas ionized by the WR/OB cluster. If this 
idea is confirmed, the H 11 region powered by the W43 
starburst cluster would cover up to 13 pc. 

Fig. 6 displays the infrared emission mapped with MSX 
at 21 jim.. It nicely coincides with the free-free emission 
of the giant H 11 region (compare with Fig. 3). Since very 
small grains (Desert, Boulanger, & Puget 1990) survive 
inside H 11 regions, the 8 to 21 ^im MSX maps may trace 
their strong near- to mid-infrared continuum. Indeed, this 
hot continuum has already been detected in a region with 
similarly intense radiation field (namely M17), using ISO- 
CAM/CVF spectra at 5 - 16 ^m and MSX data (Cesarsky 
et al. 1996; Nielbock et al. 2001). The 21 fim emission in 
Fig. 6 peaks at the location of the stellar cluster but also 
extends along the southern and northern clouds containing 
W43-MM4, MM11, MM14 and W43-MM19, MM6, MM8, 
MM32. Those ridges of mid- infrared emission may outline 
the ionization fronts also observed at centimeter and far- 
infrared wavelengths (cf. Lester et al. 1985). The relative 
positions of the 21 [im ridges to the 1.3 mm dense clouds 
perfectly agree with the 3D model of W43 (cf. Fig. 6). The 
southern 21 /im emission is displaced to the north of the 
clouds associated with W43-MM4, MM11 and MM14, sug- 
gesting that most of it is obscured by these clouds where 
we measure up to A v <~ 300 mag. In contrast, the north- 
ern 21 /jm ridge is consistent with the illumination of the 
clouds surface from the front and south. The 10 5 L Q lu- 
minosity and ~ 1 pc extent of these ridges also agree with 
the radiation of the cluster being intercepted within rea- 
sonable opening angles of 1/20 x Air sr. It is thus very 
likely that the starburst cluster, located in between these 
clouds (see e.g. Fig. 7a), is responsible for the large mid- 
infrared emission mapped in e.g. Fig. 6. Therefore, the 
two IRAS sources 18450-0200 and 18449-0158 (from the 
Point Source Catalogue) could arise at the interface be- 
tween the ionized gas and the molecular cloud rather than 
from any embedded stellar population (see also Mathis 
1990). This assumption holds when comparing the IRAS 
fluxes with those measured by MSX and the KAO which 
both have better angular resolution (20" and 50" vs. 30" 
and 60" - 120"). The 12 and 25 fim fluxes of IRAS 18450- 
0200 and IRAS 18449-0158 are indeed twice stronger than 
those of the compact sources MSX G30. 7585-0.0495 and 
MSX G30.7815-0.0219 (Egan ct al. 2001). Similarly, the 
60 and 100 /im IRAS fluxes are ten times those measured 
with the KAO (Figs, la-b in Lester et al. 1985). A near- 
infrared study will uniquely determine if these two IRAS 
sources are indeed associated with ionization fronts excited 
by the WR/OB cluster (Petr-Gotzens et al. in prep.). 

The existing data suggest that the centimeter and near- 
to far-infrared characteristics of W43 are completely dom- 
inated by the ionization and heating from the WR/OB 
cluster. We hereafter investigate if the kinematics of the 
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molecular clouds observed in Figs. 1-2 is itself shaped by 
this starburst cluster. 

4.1.2. Impact of the starburst cluster on the clouds 

To distinguish between the influence of Galactic motions 
and that of the starburst cluster, we need to carefully 
study the systemic velocity throughout the W43 cloud 
complex. Fig. 7a displays the v LSR measured for the com- 
pact cloud fragments of Table 2 as a function of the Galac- 
tic longitude offsets relative to that of the starburst cluster. 
It shows a clear velocity gradient along the Galactic plane, 
which is consistent with that mentioned in Sect. 3.1.2. The 
dashed line in Fig. 7a represents a linear fit, which gives a 
gradient of ~ 0.4 kms -1 pc -1 . Except in the improbable 
case where the W43 molecular cloud is along the line of 
sight (i < 3°), this velocity gradient is too large to be ex- 
plained by a simple model of spiral Galactic motions (e.g. 
Brand & Blitz 1993). In contrast such a gradient agrees 
with W43 being in the molecular ring, at a location where 
the velocity field is known to be complex (Liszt et al. 1993; 
Liszt 1995). 

In Fig. 7b, the v LSR is corrected for the Galactic gradient 
measured in Fig. 7a and plotted as a function of the pro- 
jected distance to the WR/OB cluster. This distance is set 
to a negative value for cloud fragments situated south of 
the stellar cluster. The dispersion of the cloud fragments 
and the velocity dispersion of individual fragments are ob- 
served to increase in the inner 4 pc from the WR/OB 
cluster, outlined by a dashed ellipse in Fig. 7b. Fur- 
thermore, the dispersion from the median velocity (zero 
level indicated as a gray line) is larger for lower density 
unbounded fragments. Indeed, the four velocity compo- 
nents detected at location W43-MM15&20 have density 
lower than 10 5 cm~ 3 and are far from virial equilibrium 
(M smn /M ti[ <~ 0.04). Besides, the high velocity compo- 
nents of W43-MM4 and W43-MM7 probably correspond 
to the lowest density clouds at the location of these dust 
fragments (cf. Sect. 3.2.3). Therefore, we find that the dis- 
persion of the fragments systemic velocities is decreasing 
with the distance to the WR/OB cluster and with the den- 
sity of the clouds. It can be interpreted as evidence that 
low- to moderate-density clouds (< 10 5 cm~ 3 , < 50 M Q ) 
situated within a 3 — 4 pc radius from the starburst clus- 
ter have been blown away. If we assume a typical age of 
10 6 yr for the WR and OB stars, the ~ 30 kms -1 dis- 
persion recorded on the line of sight of the cluster and 
supposed to correspond to the 7 pc diameter measured in 
Fig. 7b suggests a sweeping motion of ~ 3.5 kms -1 . A 
similar study around the M17 giant H n region suggests 
that a shell of cloud is radially expanding with ~ 10 km s _1 
(Rainey et al. 1987). Such velocities can be reached by the 
combined effect of stellar winds and ionization shocks orig- 
inating from the starburst cluster, since models predict a 
rocket velocity up to 5 kms -1 (Bcrtoldi 1989). 

Several studies already suggested the presence of a shell 
expanding away from the WR/OB cluster (Liszt et al. 
1993; Balser et al. 2001). We confirm this idea but pro- 
pose that the densest parts of the molecular cloud complex 
have a systemic velocity which is given by local Galactic 
motions and is not affected by the interaction with the 
WR/OB cluster. This should be particularly true for the 
vast majority of the compact cloud fragments identified 



in Sect. 3.1.1. Their detailed structure and kinematics 
can, however, be sensitive to the shock wave propagating 
ahead of the ionization front. The kinematics of the W43 
molecular cloud complex would merit a careful compari- 
son with models of "radiation-driven implosion" (Bcrtoldi 
1989; Lefloch & Lazareff 1994). Such a study will however 
need further molecular line observations and is clearly be- 
yond the scope of the present paper. 

4.2. Protocluster Candidates in the Mini- starburst W4S 

We have identified 48 cloud fragments which we think 
are good sites for harboring on-going massive star for- 
mation. Their general characteristics are investigated in 
Sect. 4.2.1 and suggest they are protoclusters, i.e. molecu- 
lar clouds that will form clusters of stars. We then discuss 
in Sect. 4.2.2 the evolutionary state of several of these pro- 
toclusters and qualify W43 as a Galactic mini-starburst in 
Sect. 4.2.3. 

4.2.1. Protocluster candidates 

In order to determine the nature of the dense, mas- 
sive, and compact cloud fragments of W43, we com- 
pare their general characteristics with those of cloud frag- 
ments observed in nearby (low- to intermediate-mass) star- 
forming regions. Those cloud fragments are frequently 
called "dense cores", according to a treshold density (> 
10 4 cm~ 3 ) and regardless of their size, mass or substruc- 
ture (e.g. Jijina, Myers, & Adams 1999). Instead, we 
prefer to use the terminology recently established for star- 
forming molecular clouds like p Ophiuchi or Orion B when 
they are surveyed with high-density and high-resolution 
tracers such as the submillimcter continuum (e.g. Motte 
& Andre 2001b). The main advantage is to suggest the 
level of subfragmentation expected for any cloud fragment 
and thus its abilility to form either one or several stars. 
Molecular clouds which are currently forming clusters of 
stars are often called protoclusters, like the p Ophiuchi 
main cloud L1688 or NGC 2068 in Orion B (e.g. Wilking 
& Lada 1983; Motte et al. 2001). The best studied case is 
the L1688 protocluster which has a mass of ~ 500 M Q 
in ~ 1 pc diameter (Wilking & Lada 1983). In those 
protoclusters, starless condensations are believed to repre- 
sent the last stage of turbulent fragmentation (~ 0.5 M Q , 
~ 0.03 pc), making them the precursors of individual stars 
(Motte et al. 1998, 2001). Dense cores are the intermediary 
cloud structures most probably forming a group of stars 
(~ 5 Mq, ~ 0.1 pc, e.g. Jijina et al. 1999). 

The W43 fragments have diameters ranging from 0.09 pc 
to 0.56 pc and masses spanning the range 20 M Q to 
3 600 M Q (cf. Tables 1 and 3). On average (~ 300 M Q , 
~ 0.25 pc), their characteristics are thus closer to those 
of protoclusters than those of low-mass dense cores or 
condensations. In that context, most of the cloud frag- 
ments observed in W43 qualify as protocluster candidates. 
Higher-resolution studies are necessary to prove our inter- 
pretation but we give below additional arguments. In- 
deed, the HCO + line characteristics of the W43 compact 
cloud fragments suggest that they are not homogeneous 
reservoirs of gas. Indeed, optically thin lines are broad 
(Av ^2 — 8 kms -1 , cf. Table 2) and often non-Gaussian 
(see e.g. Figs. 4a-e). Moreover, H 13 CO + which is tracing 
denser gas than HCO + has narrower lines which some- 
times peak at a slightly different velocity (e.g. W43-MM4, 
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MM7, MM11). Jijina et al. (1999) showed that cloud frag- 
ments associated with clusters have larger line widths than 
those forming single stars. Myers (1998) also predicted 
the existence of ~ 0.03 pc condensations embedded within 
massive and turbulent cloud fragments. Because the com- 
pact fragments of W43 are indeed massive (~ 300 M ) 
and turbulent (Av ~ 5 kms -1 ), they are likely to consist 
of many smaller-size and denser structures orbiting within 
the same cloud. This interpretation argues for the pro- 
tocluster nature of the submillimcter fragments we have 
identified. 

With present observations only, we cannot draw definite 
conclusions on the gravitational status of W43 protoclus- 
ters (see Sect. 3.2.3). The small M Bmm /M vir ratios observed 
for some of the most massive protocluster candidates sug- 
gest the majority of them could be far from virial equilib- 
rium. Besides, the mean ellipticity of the W43 protoclus- 
ters (axis ratio 1:2) resembles that observed for low-mass 
dense cores and is definitively larger than that of gravita- 
tionally bound starless condensations (Motte et al. 2001). 
It thus agrees with the idea that W43 protoclusters are 
generally not in a gravitation-dominant regime but prob- 
ably consist of material governed by turbulence. However, 
protocluster candidates are embedded within some hot- 
ter lower-density gas, so that the external pressure could 
substantially help confine the protocluster gas and maybe 
prevent it from dispersing. Note that the structures we de- 
fine here as protocluster candidates do not need to evolve 
independently from their surrounding gas, in which one 
expects low- to intermediate-mass stars to form. 

Strikingly, these protocluster candidates are very dense 
(mean density of n H2 ~ 10 6 cm -3 , see Table 3) compared 
to their low-mass counterparts ( ~ 10 4 cm~ 3 , e.g. Wilking 
& Lada 1983). In fact they are as dense as the low-mass 
starless condensations, which are believed to be the direct 
progenitors of single (or binary) stars (Motte et al. 1998). 
Such a dense cloud medium is thus regularly observed in 
low-mass star-forming clusters, but not over ~ 0.25 pc 
scale, which is ten times larger than the typical size of 
starless condensations. Furthermore, we stress that the 
density of smaller-scale fragments within those protoclus- 
ters will surely be much higher. Similar densities, sizes, 
masses and line widths have already been recorded for pro- 
toclusters harboring UCH us and HMPOs (Hunter et al. 
2000; Plume et al. 1999; Beuther et al. 2002). Such cloud 
characteristics may then enhance star formation and allow 
the formation of massive stars. 

The unbiased survey we made in the W43 main complex 
should be sensitive to protoclusters either being pre-stellar 
or containing massive protostars or UCH n regions. As all 
sources are at the same distance from the Sun and in the 
same cloud complex, the resulting sample is far more ho- 
mogeneous than in any other previous study (e.g. Hunter 
et al. 2000; Sridharan et al. 2002). No doubt that follow- 
ups of such a complete and comprehensive sample will help 
understand the formation mechanism of massive stars. We 
present below our first attempt to determine the evolution- 
ary state of W43 protoclusters candidates. 

4.2.2. Evolutionary state of W43 protoclusters 



It is questionable to try to determine the evolutionary 
state of protoclusters since, by definition, they contain sev- 
eral YSOs. However when a high-mass star is forming, it 
probably shapes the structure, kinematics, and chemistry 
of the whole protocluster via strong heating, ionization 
and outflow. It is thus tempting to constrain the evo- 
lution of protoclusters using the evolutionary state of its 
main stellar component. We discuss below the case of 6 
protocluster candidates according to their characteristics 
and the stellar activity signposts presented in Sect. 3.2. 

W43-MM1, MM2 and MM11 are methanol, water (with 
the exception of MM2) and hydroxyl maser sources with- 
out any infrared or centimeter detection (cf. Figs, lb, 3 
and 6). W43-MM1 and MM2 are also massive (M Bmm > 
1000 Mq) luminous (L bol > 10 4 L@) protoclusters, which 
are likely to be gravitationally bound (M smm /M vir > 3). 
The gravitational status of W43-MM11 is less clear, since 
we measure M smm /M vir ~ 0.35 (see Sect. 3.2.3). These 
protoclusters are globally cold (T dust ~ 20 K, e.g. Fig. 5a), 
but recent observations argue for the presence of embed- 
ded "hot cores" with T > 200 K in W43-MM1 and MM2 
(Motte et al. in prep.). The submillimeter to bolometric 
luminosity ratio measured in Sect. 3.2.1 (L x > 350 Mm /£ b oi ~ 
1.5 — 3%) is similar to that of low-mass class protostars 
which are YSOs in their main accretion phase (Andre, 
Ward-Thompson, & Barsony 2000). All those rare char- 
acteristics make these three protoclusters excellent can- 
didates for containing massive protostars, i.e. massive 
YSOs which have not yet assembled the bulk of their fi- 
nal stellar mass and neither have developed H n regions. 
As they are not coincident with any near- to mid-infrared 
sources, they likely represent even younger phases than 
the IRAS selected HMPOs (Brand et al. 2001; Sridha- 
ran et al. 2002). This idea is confirmed by their high 
mass to luminosity ratio, sometimes used as an evolution- 
ary indicator. W43-MM1 and MM2 have higher ratios 
(M smm /i bol ~ 0.1 - 0.2 Mq Lq- 1 ) than those measured 6 
for HMPOs (~ 0.01 M & L Q ~\ Sridharan et al. 2002) and 
UCH us envelopes (~ 0.006 M QJ L© _1 , cf. Hunter et al. 
2000). 

The W43-MM3 protocluster looks very similar to W43- 
MM1 and MM2 in terms of mass, luminosity, average 
temperature and other mass or luminosity ratios (see Ta- 
ble 3 and Sect. 3.2.1). However, W43-MM3 is not a maser 
source and coincides with a bright (S 3 . 5 cm ~ 1 Jy) ra- 
dio continuum source which is partially optically thick be- 
tween 21 and 3.5 cm (see Fig. 3 and Sect. 3.2.4). This 
compact free-free emission source is barely resolved, with 
a ~ 0.13 pc deconvolved diameter smaller than the 0.25 pc 
size measured at 1.3 mm. Despite its high luminosity 
(< 10 4 L Q ), W43-MM3 was not detected by IRAS, maybe 
due to confusion from the closeby giant H n region. The 
inspection of Fig. 6 suggests that W43-MM3 is loosely as- 
sociated with MSX G30. 7196-0.0854, but this source is 
three times more extended than the 1.3 mm protoclus- 
ter. Moreover, its spectral index measured from 8 /im to 
21 fxm indicates that the MSX source is not an embedded 
star/cluster (Lada 1987, see also Fig. 5b), but probably 
lies by chance in front of the W43-MM3 protocluster. We 
thus propose that W43-MM3 harbors an ultra-compact 



6 In order to ensure consistency with our derivation, the mass values of Beuther et al. (2002) and Hunter et al. (2000) have been recalculated 
using Ki.3 mm = 0.01 cm 2 g _1 . 
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H II region, not developed enough to efficiently heat its 
~ 1 000 M Q of dust and gas. This source illustrates the 
difficulty and inadequacy of the determination of the mas- 
sive YSOs evolutionary state from the characteristics of 
their hosting protocluster and vice versa. 

W43-MM4 is a massive protocluster (M Bmm ~ 500 M ), 
which is probably gravitationally bound (M Bmm /M vir ~ 
1.4— 0.7, cf. Sect. 3.2.3). Its study at infrared and centime- 
ter wavelengths is unfortunately impaired by the strong 
emission of the southern ionizing front (see Sect. 4.1.1). 
However in the 21 cm and 3.5 cm surveys of Zoonematk- 
ermani et al. (1990) and Becker ct al. (1994), W43-MM4 
coincides with some compact (~ 0.09 pc) free-free emis- 
sion, which is strong (S 6 cm ~ 300 mJy) and partially opti- 
cally thick, i.e. typical of UCH n regions. High-resolution 
surveys at centimeter and mid-infrared wavelengths are, 
however, necessary to firmly identify W43-MM4 as a pro- 
tocluster hosting an UCH II region. 

The W43-MM13 protocluster has a mass of ~ 200 M 
and coincides with IRAS 18450-0205 (or MSX G30.6882- 
0.0726) and some free-free emission peaks detected at 
6 cm and 21 cm (see Figs. 6 and 3). The IRAS source 
is satisfying the criteria established by Wood & Church- 
well (1989) for selecting UCH n regions by their mid- to 
far-infrared colors (see HIRES flux in Sect. 3.2.4). How- 
ever, both the MSX and centimeter sources are extended 
{FWHM ~ 1.9 pc at 21 cm and 0.9 pc at 21 urn), rul- 
ing out the interpretation of W43-MM13 containing an 
UCH II. Since the spectral index measured from 8 /im 
to 60 /jm is characteristic of a deeply embedded stellar 
source (Lada 1987), we suggest that a well developed 
H II region is observed behind a foreground cloud. The 
visible to mid-infrared radiation of such an H n region 
(L bol ~ 10 5 Lq) could easily be extinguished and red- 
dened by the A v ~ 80 mag we measure for the cloud 
complex containing W43-MM13. Another alternative is 
that the IRAS/MSX and VLA 3 cm/21 cm sources mark 
a third ionizing front excited by the WR/OB cluster (see 
Sect. 4.1.1). 

We do not have enough information to determine the 
nature of the remaining protoclustcrs of Table 3. Follow- 
up observations, searching for "hot core" and outflow sig- 
natures, as well as ultracompact free-free emission are in 
progress. Mid- to far-infrared observations with the forth- 
coming SIRTF and Herschel satellites will also help in de- 
termining their evolutionary state. 

4.2.3. The mini-starburst region W^3 

As illustrated above, the main molecular complex of 
W43 looks very efficient in forming massive stars. In- 
deed, at this early stage, we have already identified five 
protoclusters containing at least one massive protostar or 
UCH II (Sect. 4.2.2). The - 15 protoclusters in Table 3, 
which have high density (e.g. n H2 ~ 5 x 10 5 — 9 x 10 6 cm~ 3 ) 
and high mass (e.g. M Bmm ~ 100 — 3 600 Mq) are also 
good candidates for being sites of on-going or future mas- 
sive star formation. Interestingly, our results preclude that 
star formation in the W43 main molecular complex is con- 
tinuous with an efficiency typical of giant molecular clouds 
(SFE~ 1% over 10 7 yr, Silk 1997). Indeed with a Salpeter 
(1955) initial mass function, one expects to discover less 
than one massive YSO in the HMPO or UCH n phases 



(lasting for 10 4 — 10 yr) when mapping the W43 main 
complex. Therefore, we speculate that star formation is 
synchronized within the W43 main molecular complex. If 
the ~ 15 protoclusters identified above are indeed cur- 
rently forming or on the verge of forming massive stars, the 
10 6 Mq giant molecular cloud is likely to disperse within 
~ 10 6 yr, which is the lifetime of a massive star. Dur- 
ing this short period, stars should form with an efficiency 
that could be as high as SFE ~ 25%, or equivalently a 
star formation rate as high as ~ 0.25 MQyr -1 . Such a 
star formation efficiency is one order of magnitude larger 
than that of normal giant molecular clouds. It is similar 
to that found for the p Ophiuchi dense cores, whose size is 
one hundred times smaller than the W43 main molecular 
complex (SFE~ 31% over 10 6 yr, Bontemps et al. 2001). 
In 10 6 years from now, if W43 becomes a complex of OB 
clusters encompassed in the same (14 pc) 3 volume, this 
efficiency would correspond to a stellar density of ~ 100 
stars/pc 3 (assuming a Salpeter IMF). A much larger den- 
sity has been recorded in the Arches cluster and the core of 
the NGC 3603 starburst cluster (10 5 stars/pc 3 , e.g. Tapia 
et al. 2001; Figer et al. 1999) but over a 10 4 - 10 6 times 
smaller volume. Although this result clearly needs to be 
confirmed, it suggests that the W43 molecular complex 
is experiencing a sudden burst of star formation. There- 
fore, W43, which is known to contain a starburst cluster 
of main sequence stars, probably also harbors a molecular 
region still actively forming massive stars. It thus very 
likely qualifies as a Galactic mini-starburst region, a term 
we coin to define a region whose cloud and stellar content 
could be used as a template for those in starburst galaxies. 

Since W43 is undergoing a mini-starburst, it offers the 
unique opportunity to determine the physical processes re- 
sponsible for enhancing star formation in molecular cloud 
complexes. Among our results we show that, despite the 
extremely high radiation emitted by the WR/OB cluster, 
the densest parts of the W43 molecular clouds remain cool. 
Such low temperatures (T dust ~ 20 K) are unusual for sites 
of massive star formation where T dust > 40 K is measured 
(Hunter et al. 2000; Sridharan et al. 2002). Low tem- 
peratures coupled with high densities, could be one key 
element of the recipe for local mini-starbursts. Another 
obvious element may be the impact of the revealed clus- 
ter of WR/OB stars onto the cloud. Further studies are 
requested before one can make any breakthrough on the 
origin of such mini-starburst in our Galaxy. 

5. SUMMARY AND CONCLUSION 

W43 is a molecular and H II complex whose ionized gas 
has been far better studied than its molecular component. 
We have imaged the W43 main cloud at submillimeter 
wavelengths in the 1.3 mm and 350 fim continuum and 
HCO + (3-2) line emission. In addition, we have obtained 
deep HCO+(3-2) and H 13 CO+(3-2) spectra at selected lo- 
cations. Our main findings can be summarized as follows: 

1. A multircsolution analysis on our submillimeter 
continuum maps identifies ~ 50 compact fragments. 
Our 350 ^.m continuum observations along with the 
3.5 cm image of Balser et al. (2001) show that their 
1.3 mm emission is largely thermal emission from 
cool dust. 
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2. These bona- fide cloud fragments have diameters 
varying from 0.09 pc to 0.56 pc and masses span- 
ning the range 20 M© to 3 600 M . Their large 
size and turbulent line width (Aw ~ 5 kms -1 ) sug- 
gest they are protoclusters, i.e. clouds that con- 
tain many smaller-size and denser structures and 
will form star clusters. Those protocluster candi- 
dates have large mean densities (n H2 ~ 10 6 cm~ 3 ) 
reminiscent of the direct progenitors of individual 
low-mass stars. The W43 protocluster candidates 
thus constitute an excellent sample for studies of 
the earliest stages of massive star formation. 

3. The present, unbiased survey is sensitive to proto- 
clusters either being pre-stellar or containing mas- 
sive protostars or UCH II regions. Five of the 
W43 protoclusters are confirmed to contain massive 
YSOs in their HMPO or UCH II phase. Follow-up 
observations are needed to determine the evolution- 
ary state of the remaining protocluster candidates. 

4. In W43, the centimeter and infrared emission from 
the giant H II region dominates, preventing the de- 
tection of more compact sources that could be as- 
sociated with UCH us. Notably, at least two of the 
three IRAS point sources of the W43 main cloud 
are not stellar in nature; they are instead associ- 
ated with ionization fronts possibly all excited by 
the closeby WR/OB cluster. 

5. While the low-density clouds surrounding the star- 



burst cluster may have been blown away, the dens- 
est parts of the W43 molecular clouds seem to re- 
main at the systemic velocity given by Galactic mo- 
tions. A more precise kinematic study is necessary 
to determine if the ionization front associated with 
the giant H n region is compressing the molecular 
clouds and triggering star formation. 

6. W43 is the site of at least two remarkably efficient 
episodes of massive star formation. Indeed, it is 
known to harbor a starburst cluster containing sev- 
eral WR and OB stars (~ 10 6 L Q ). We show here 
that the molecular complex is currently undergoing 
a second mini-starburst with a star formation effi- 
ciency of ~ 25%/10 6 yr and possibly a final stellar 
density of ~ 100 stars/pc 3 over (14 pc) 3 . Learn- 
ing about the global characteristics of this Galactic 
mini-starburst region should help constraining the 
physical processes at work in the distant starburst 
galaxies. 
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a (2000) a (2000) 

Fig. 1. — The W43 main star-forming complex mapped at 1.3 mm with MAMBO (a) and 350 fim with SHARC (b). Contour levels are 
logarithmic, spaced by factors 2. Levels and rms noise are: (a) 60, 120, 240, 480, 960 mjybeam -1 , and 2, 4 Jybeam" 1 , la ~ 12 mjy/ll"- 
beam; (b) 12, 24, 48, 96 and 200 Jybeam" 1 , lcr ~ 3 Jy/11.5"-bcam. In (a), the G aussian dust fragments extracted from the 1.3 mm map 
are indicated by ellipses, fragments which are relevant for the discussion arc also labelled. In (b), the WR/OB association is indicated by a 
star symbol. OH, H2O and CH3OH maser sources are also shown as open squares, filled triangles, and three branches crosses, respectively. 
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Compact fragments detected at 1.3 mm and 350 /j,m 
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a Deconvolved FWHM size derived from a 2D-Gaussian fit to the 1.3 mm map (after background subtraction, cf. 
Sect. 3.1.1). An upper limit FWHM of 0.09 pc was assumed for unresolved fragments. 

b Integrated flux measured at 350 /xm (after background subtraction). Upper limits of 2 x 3a are given when the 
1.3 mm fragment is undetected at 350 fim. The factor 2 accounts for the average SV.^mm/^i 1 ^ mm ratio of these compact 
but resolved sources. 

c Free-free contamination of the 1.3 mm fluxes estimated from the 3.5 cm map of Balser et al. (2001) or the 21 cm 
map of Liszt et al. (1993) when indicated by a star maker. 
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Fig. 2. — The W43 main star-forming complex mapped in HCO + (3-2). In (a), the intensity map is integrated over the 84 — 104 kms -1 
velocity range. Contour levels are 9 to 36 by 9 Kkms -1 and rms noise is la ~ 3 Kkms -1 . In (b)-(f), the channel maps are integrated 
over 4 km s _1 channels. Contour levels arc 0.7 to 4.2 by 0.7 Kkms" 1 and rms noise is la ~ 0.23 Kkms 1 . The positions of dust fragments 
extracted from Fig. la are marked by crosses. The Galactic plane at the WR/OB cluster latitude is displayed in (e). 
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18 h 47 m 40 s 
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Fig. 3. — The spectral index map of W43 measured between 350 (im and 1.3 mm (gray scale) compared with the free-free emission mapped 
at 21 cm with the VLA (contours, cf. Liszt ct al. 1993). The submillimctcr spectral index varies from 2 (black level) to 4 (white 

level) with a step of 0.5. Contour levels are 10% to 100% by 10% with a maximum centimeter flux of ~ 0.8 Jy/12.5"-beam. The WR/OB 
association is indicated by the star symbol, the ionization front caused by its UV radiation is outlined by a thick dashed white ellipse. 
Small-diameter sources identified at 21 cm, 6 cm and 3.5 cm wavelengths are marked by filled squares, circles, and triangles, respectively. 
Selected submillimctcr fragments are indicated with crosses. 
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Fig. 4. — HCO+(3-2) (thin line histogram) and H 13 CO + (3-2) (thick line histogram) spectra smoothed to a resolution of 0.6 kms 
(a)-(b) and (d)-(e) the sources are located north (resp. south) of the WR/OB cluster which lies at location W43-MM15&20. 
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Fig. 5. — Spectral energy distributions of the cloud fragments W43-MM1 (a) and W43-MM3 (b) compared with gaybody models with 
varying dust opacity indices and temperatures. The absolute uncertainty of MSX and KAO fluxes is set to 30%, that of MAMBO and SHARC 
fluxes is 20%. 
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Fig. 6. — The continuum emission of W43 at 21 (im (gray scale) compared with that at 1.3 mm (contours). The 21 fim map was obtained 
by the MSX satellite with a 20" aperture. Levels are 0.5, 1, 3, 6, and 12 xlO 9 Jysr~ x . The 1.3 mm map is shown in Fig. la with the 
same contours, except the first one. IRAS/MSX sources and the WR/OB association are indicated by star symbols; selected submillimctcr 
fragments are labelled. 
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Fig. 7. — Systemic velocity of selected submillimctcr fragments as a function of the Galactic longitude offsets (a) and projected distance 
(b) from the WR/OB cluster. The dashed line in (a) marks the least squares fit gradient which is subtracted from the vlsr plotted in (b). 
The location of the WR/OB cluster is indicated with a star marker. Error bars correspond to the line width of the velocity components. Note 
that in (b), the dispersion to the zero line (gray thick line) and the line widths are larger within the dashed ellipse. 
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Table 2 





HCO + 


LINE SURVEY AT 


SELECTED 


LOCATIONS 






rragmcnt 


Molecular 


Component 


^LSR 


± MB 


J 1 MB av 


A n i 


name 


transition 


number 


|kms J 


[KJ 


[K kms J 


[kms J 


W43-MM1 


H 13 CO+(3-2) 


1 


98.8 


1.30 ±0.10 


8.2 


5.9 


W43-MM2 


H 13 CO+(3-2) 


1 


90.8 


1.04 ±0.09 


4.6 


4.2 


W43-MM3 


H 13 CO+(3-2) 


1 


93.5 


0.82 ±0.06 


3.8 


4.4 


W43-MM4 


HCO+(3-2) 


1 


91.8 


3.24 ±0.14 


17.5 


5.1 




2 


97.6 


1.83 ±0.14 


6.0 


3.1 




H 13 CO+(3-2) 


1 


92.0 


1.55 ±0.13 


4.4 


2.7 






2 


97.3 


0.20 ±0.13 


0.7 


3.2 


W43-MM6&8 


HCO+(3-2) 


1 


94.4 


2.37 ±0.14 


20.9 


8.3 




H 13 CO+(3-2) 


1 


94.3 


0.35 ±0.06 


2.4 


6.3 


W43-MM7 


HCO+(3-2) 


1 


90.6 


0.83 ±0.14 


3.9 


4.4 




2 


97.4 


1.72 ±0.14 


6.7 


3.7 




H 13 CO+(3-2) 


1 


91.9 


0.19 ±0.06 


0.7 


3.7 






2 


96.0 


0.31 ±0.06 


0.7 


2.2 


W43-MM9 


H 13 CO+(3-2) 


1 


96.0 


0.61 ±0.11 


3.6 


5.5 


W43-MM10 


H 13 CO+(3-2) 


1 


92.9 


0.59 ±0.11 


2.3 


3.7 


W43-MM11 


HCO+(3-2) 


1 


91.5 


2.02 ±0.10 


14.4 


6.7 




H 13 CO+(3-2) 


1 


91.8 


0.80 ±0.13 


3.2 


3.8 


W43-MM14 


HCO+(3-2) 


1 


89.9 


1.01 ±0.08 


6.6 


6.1 


W43-MM15&20 


HCO+(3-2) 


1 


81.0 


0.48 ±0.05 


2.4 


4.8 




2 


87.0 


0.23 ±0.05 


0.9 


3.5 






3 


94.4 


0.58 ±0.05 


4.2 


6.8 






4 


107.7 


0.10 ±0.05 


0.7 


6.0 


W43-MM25 


HCO+(3-2) 


1 


92.3 


0.67 ±0.09 


3.5 


4.9 
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Table 3 

Properties of submillimeter compact fragments 



Fragment 


T dust 






> a 




nb xM viT b 


name 


[K] 


Wq] 


[cm" 


3 ] 


[Mq] 


"VT T 7 A i~> TV TH ,T -| 

W43-MM1 


20 


3590 


Q C 
O.O 


X 


10 b 


1040 


TTT^f* "TV /T "IV 1~ i~\ 

W43-MM2 


20 


1620 


A ft 
4.0 


X 


10 6 


490 


"VT T 7 A C\ TV /T TV TO 

W43-MM3 


20 


960 


z.U 


X 


10 6 


610 


IT- t , \ ~\ IT~\ ,T A 

W43-MM4 


20 


520 




X 


10 6 


2 x 370 


W43-MM5 


20 


110 


O.O 


X 


10 6 




W43-MM6 


20 


500 


o n 
z.U 


X 


10 5 


2930* 


it- 1 . "\ TV .1" TV T" 

W43-MM7 


20 


870 


1 7 
1. / 


X 


10 5 


2 x 840 


it- 1 . \ TV .1" TV .1" 

W43-MM8 


20 


390 


1 
l.Z 


X 


10 5 


2930* 


W43-MM9 


20 


290 


0.0 


X 


10 5 


920 


"VT T 7 < f» TV /I" TV /I" 1 1""1 

W43-MM10 


20 


230 


1 s 
l.o 


X 


10 6 


270 


W43-MM11 


20 


190 


1.0 


X 


10 6 


530 


"Y T 7 * n TV /TTV /I" 1 i~\ 

W43-MM12 


20 


90 


1 7 
1. 1 


X 


10 6 




W43-MM13 


20 


200 


9 K 


X 


10 5 




TIT in TV /I" TV /I" -| 1 

W43-MM14 


20 


230 


Q Q 
O.O 


X 


10 5 


1260 


W43-MM15 


30 


200 


1 1 
1 . 1 


X 


10 5 


4 x 1380 


"T T 7 *f n 1\ Tl T-| 

W43-MM16 


20 


260 


A 9 
4.Z 


X 


10 5 




W43-MM17 


20 


310 


9 K 
Z.O 


X 


10 5 




W43-MM18 


20 


180 


7 £ 
I .0 


X 


10 5 




T T 7 < f"» TV /T TV IT 1 l"l 

W43-MM19 


20 


180 


i n 

l.U 


X 


10 5 




W43-MM20 


30 


50 


1 

i.y 


X 


10 5 


4 x 1380* 


"T T 7 ^ n TV /T TV /T i~\ -I 

W43-MM21 


20 


210 


9 n 


X 


10 5 




W43-MM22 


20 


320 


i q 

l.O 


X 


10 5 




Til in TV /TTV /I" n 

W43-MM23 


20 


140 


/I Q 

4.y 


X 


10 5 




^T7^^ x ^ t\ /r t\ /r i~\ < 

W43-MM24 


20 


370 


y.i 


X 


10 4 




ITT i n TV ,t "ft ,r r- 

W43-MM25 


30 


50 


O.O 


X 


10 5 


490 


W43-MM26 


20 


230 


i /i 

1.4 


X 


10 5 




"T T 7 < n TV. /TT\ /T i~\ 

W43-MM27 


20 


110 


/I R 
4.o 


X 


10 5 




^T7^n t\ /r t\ /r i~\ r-i 

W43-MM28 


20 


110 


o.y 


X 


10 5 




TT7^n TV r TV T iT\ l"l 

W43-MM29 


20 


90 


y .o 


X 


10 5 




W43-MM30 


20 


300 


O.O 


X 


10 4 




tl7*n tv r tv r n 1 

W43-MM31 


20 


130 


9 ft 
z.O 


X 


10 5 




it- i . \ i ,ti\ fiiii 

W43-MM33 


20 


110 


0.1 


X 


10 5 




W43-MM34 


20 


70 


4.U 


X 


10 5 




W43-MM35 


30 


110 


£ 9 
O.z 


X 


10 4 




W43-MM36 


30 


40 


o.y 


X 


10 5 




W43-MM37 


30 


40 


4.0 


X 


10 5 




it- i . \ -iv /r-\ ,r n 

W43-MM38 


30 


140 


9 Q 
Z.O 


X 


10 4 




^T7^n tv r tv r n t~\ 

W43-MM39 


20 


120 


1.8 


X 


10 5 




W43-MM40 


30 


60 


8.5 


X 


10 4 




~T~t T A i~> J ITU f 1 1 

W43-MM41 


30 


80 


3.4 


X 


10 4 




W43-MM42 


20 


60 


3.9 


X 


10 5 




W43-MM43 


30 


50 


1.7 


X 


10 5 




W43-MM44 


20 


100 


1.0 


X 


10 5 




W43-MM45 


30 


20 


1.0 


X 


10 6 




W43-MM47 


30 


50 


8.9 


X 


10 4 




W43-MM48 


20 


60 


3.0 


X 


10 5 




W43-MM49 


20 


50 


2.2 


X 


10 6 




W43-MM51 


20 


40 


1.8 


X 


10 6 





a Mean density derived from Col. 3, and Col. 5 of Ta- 
ble 1: < n H2 >= M smm / [§tt x (FWHM/2) 3 ]. 

b Virial mass of an average velocity component multi- 
plied by the number of components. Virial masses that 
correspond to two independent submillimeter fragments 
are indicated by a star marker. 



